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Abstract 

We describe a general multiplier method to obtain boundary stabilization of the wave equation by 
means of a (linear or quasi-linear) Neumann feedback. This also enables us to get Dirichlet boundary 
control of the wave equation. This method leads to new geometrical cases concerning the "active" 
part of the boundary where the feedback (or control) is applied. 

Due to mixed boundary conditions, the Neumann feedback case generate singularities. Under a simple 
geometrical condition concerning the orientation of the boundary, we obtain a stabilization result in 
linear or quasi-linear cases. 

Introduction 

In this paper we are concerned with control and stabilization of the wave equation in a multi-dimensional 
body fid" 

Stabilization is obtained using a feedback law given by some part of the boundary of the spacial 
domain and some function defined on this part. The problem can be written as follows 

u" - Au = in f2 x R* + , 
u = onc)l] D xR;, 
d y u = F ondn N xW + , 
u(0) = Mo m i 

u'(0) = u\ in Q , 

where we denote by u', u" , Au and d u u the first time-derivative of u, the second time-derivative of the 
scalar function u, the standard Laplacian of u and the normal outward derivative of u on <9fi, respectively; 
(dQo, OCIn) is a partition ofdfl and F is the feedback function which may depend on state (u,u'), position 
x and time t. 

Our purpose here is to choose the feedback law, that is to say the feedback function F and the "active" 
part of the boundary, DQn, so that for every initial data, the energy function 

E(t) = \ f (\u>(t)\ 2 + \Vu(t)\ 2 )dx, 

is decreasing with respect to time t, and vanishes as t — ► oo. 
Formally, we can write the time-derivative of E as follows 



E'{t) = [ Fu'da 
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and a sufficient condition for E to be non- increasing would be: Fu' < on 

Thanks to the multiplier method introduced by L.F. Ho [12] in the framework of Hilbert Uniqueness 
Method |13j . it can be shown that the energy function is uniformly decreasing as time t tends to oo by 
choosing x i— > m(x) = x — xo, where Xo is some given point in R™ and 

dn N = {xeffl/ m(x).i/(x) > } , F = -m.vu' , 

where v is the normal unit vector pointing outward of Q. 

This method has been performed by many authors, see for instance Komornik and Zuazua [llj . Komornik 
[TO] and the references therein. Here we extend the above result for rotated multipliers defined in [16] 
and we follow the analysis of singularities initiated by Grisvard [5] and extended by Bey, Loheac 
and Moussaoui @]. This last work leads to results in case of higher dimensional domains with a non- 
empty boundary interface T = dfl n H 8Qd under an additional geometrical assumption concerning the 
orientation of the boundary. 

Concerning the control problem, our goal is to find v such that the solution of 

u"-Au = infix(0,r), 

u = ondn D x (0,T), 

u = v on 90jv x (0, T), 

u(0) = u in Q, 

u (0) = u\ in Q, 

reaches an equilibrium at t = T. 

We here follow |12j : in this work, Ho used the multiplier technique. His main purpose was to prove an 
inverse inequality for the linear wave equation implying its exact controllability. He introduced the so- 
called exit condition: the control region must contain a subset of the boundary where the scalar product 
between the outward normal and the vector pointing from some origin towards the normal is positive. 
By varying the origin, a family of boundary controls satisfying the condition is obtained. 
In the last decades, micro-local techniques and geometric optics analysis allowed to find geometrical 
characterization of control and minimal control time in the exact controllability of waves. This condition 
has been introduced in [3] under the name of Geometric Control Condition ( GCC). It generalized the 
previous exit condition. 

There is a certain balance: with GCC, control time is optimal but the observability constant is not 
explicit. With exit condition, time is not optimal but observability constants can be explicit, which is 
very useful in theoretical and numerical estimations. 

In this paper we extend the family of multipliers recently introduced by Osses [16j . 

1 Notations and main results 

Let f2 be a bounded open connected set of M. n (n > 2) such that 

dCl is C 2 in the sense of Necas |15j . (1) 

In the sequel, we denote by I the n x n identity matrix and by A s the symmetric part of a matrix A. Let 
m be a W 1,00 (f2) vector-field such that 

cssinf (div(m)) > esssup(div(m) — 2A m ) (2) 

n n 

where div is the usual divergence operator and A m (x) is the smallest eigenvalue of the real symmetric 
matrix Vm(x) s . Using Sobolev embedding, one may also assume that m £ C(f2). 

Remark 1 The set of all W 1,00 (fi) vector-fields such that ^ holds is an open cone. If m belongs to 
this set, we denote 

c(m) = -(cssinf (div(m)) — ess sup (div (m) — 2A m ) 
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Examples 

• An affinc example is given by 

m(x) = (A 1 +A 2 )(x-x ), 
where Ai is a definite positive matrix, A2 a skew-symmetric matrix and xo any point in R n . 

• A non linear example is 

m(x) = {dl + A)(x- x ) + F(x) , 

where d > 0, A is a skew-symmetric matrix, Xo any point in R n and F is a W 1,00 (f2) vector field 
such that 

ess sup ||VF S || < - , 

n n 

where || • || stands for the usual 2-norm of matrices. 

We consider a partition (9J7jv, dilo) of dft such that 

T = dilu n 9f2jv is a C 3 -manifold of dimension n — 2 such that m.u = on T, 
3 oj neighborhood of T such that dil Dlu is a C 3 -manifold of dimension n — 1, (3) 
H 11 ^ 1 (OHd) > (H 71 ^ 1 is the (n — l)-dimensional Hausdorff measure) . 

Furthermore, we assume 

dn N C {xe 5fi/m(x).i/(x) > 0}, dQ D c {x e 90 / m(x).i/(x) < 0} . (4) 

This assumption clearly implies: m.u = on T. 

Boundary stabilization 

Let g : R — > R be a measurable function such that 

g is non-decreasing and 3k + > : \g(s)\ < k + \s\ a.e. . (5) 
Let us now consider the following wave problem, 

in Q x m 



where initial data satisfy 



+ ' 

u = ondfl D xR* + , 
d v u = rn.v g{u') on x W + , 
u(0) = «o in , 

u (0) = ui in SI , 



K!n)e4(Sl)xL 2 (Sl) 



with H^(fi) = {« e H^fi) : v = on 

Problem (5) is well-posed in this space. Indeed, following Komornik [TU] , we define the non- linear operator 
W on H;£,(fi) x L 2 (f7) by 

W(u,w) = (-«, - Ait) , 

fl(W) = G Hj,(l!) x H|,(ll) / An e L 2 (!l) and d v u = -m.v g(v) on dfl N } , 



so that (S) can be written as follows, 



(u,v)' + W(u,v) = 0, 
(u,v)(0) = (u ,ui) . 



It is classical that W is a maximal-monotone operator on 11^(0) x L 2 (f2) and that D(W) is dense in 
H^,(f2) x L 2 (S1) for the usual norm. Following Brczis pQ, we can deduce that for any initial data (uq, vq) in 
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D(W) there is a unique strong solution (it, u) such that u G W 1 ' O0 (R; H^(fi)) and Am G L°°(R + ; L 2 (ft)). 
Moreover, for two initial data, the corresponding solutions satisfy 

Vi > 0, IKu 1 ^),^ 1 ^)) - (u 2 (t),« 2 (t))|| H | 3 (n)xL=(n) < IIOo^'o) - ("o,«o)llHi(n)xL=(n) • 

Using the density of Z?(W), one can extend the map 

D(W) — > Hij(n) x L 2 (ft) 
(u ,w ) 1 — ► (u(t),v(t)) 

to a strongly continuous semi-group of contractions (S(t)) t >o and define for (uq,Ui) G H^(il) x L 2 (£l) 
the weak solution (u(i),t/(t)) = S(t){u ,ui) with the regularity u G C(R+; H^(fi)) n C 1 (R+; L 2 (ft)). We 
hence define the energy function of solutions by 

E (°) = I I (Kl 2 + |Vw | 2 )rfx and £(t) = - / (\u'(t)\ 2 + \Vu(t)\ 2 )dx, for t > . 

In order to get stabilization results, we need further assumptions concerning the feedback function g 

3p>l,3fc_>0: \g(s)\ > min{|s|, |s| p } , a.e. , (6) 
and the additional geometric assumption 

m.r < onT, (7) 

where r(x) is the normal unit vector pointing outward of at a point x G T when considering <9f2jv 
as a sub- manifold of dft. 

Remark 2 It is not necessary to assume that 

H n -\{x G dCl N I m(x).i/(x) > 0}) > 

to get stabilization. In fact, our choices of m imply such properties (see examples in Section 5) whether 
the energy tends to zero. 

A main tool in this work is Rcllich type relations [IT]. They lead to results of controllability and 
stabilization for the wave problem (see [TT] and [12j). When the interface V is not empty, the key- 
problem is to show the existence of a decomposition of the solution in a regular and a singular parts 
(see [7j |9]) in any dimension. The first results towards this direction are due to Moussaoui [14] . and 
Bey-Loheac-Moussaoui [J . 

In this new case, our goal is to generalize those Rellich relations. This will lead us to get a stabilization 
result about (S) under (U|), {?])■ As well as in [TO], we shall prove here two results of uniform boundary 
stabilization. 

Exponential boundary stabilization 

We here consider the case when p = 1 in j0|). This is satisfied when g is linear, 

3a > : V.s G R, g(s) = as. 
In this case, the energy function is exponentially decreasing. 

Theorem 1 Assume that conditions (JTJ) 7 ([2"]). @ and © hold and that the feedback function g satisfies 
© and © with p = l. 

Then under the further geometrical assumption ([?]). there exist C > and T > such that for every 
initial data in H^fi) x L 2 (f2), the energy of the solution u of (S) satisfies 

Vt>T, E(t) <E(0) exp(l- -) . 
The above constants C and T do not depend on initial data. 
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Rational boundary stabilization 

Wc here consider the case p > 1 and we get rational boundary stabilization. 

Theorem 2 Assume that conditions ([T]) 7 ([2"]). @ and Q /io/d and £/ia£ the feedback function g satisfies 
(0) and © ™£/i p > 1 . 

Then under the further geometrical assumption ([7]) ; there exist C > and T > such that for every 
initial data in H^(f2) x L 2 (f2), the energy of the solution u of (S) satisfies 

Vt>T, E{t) <Ct 2/(p ~ 1} . 
where C depends on the initial energy E(0). 

Remark 3 Taking advantage of the work by Banasiak and Roach JBj who generalized Grisvard's results 
in the piecewise regular case, we will see that Theorems^ and\^ remain true in the bi- dimensional 
case when assumption ([1]) is replaced by following one, 



dfl is a curvilinear polygon of class C 2 , 

each component of dil \ T is a C 2 -manifold of dimension 1 , 



and when condition ([7]) is replaced by 

Vx G T , < tu x < 7r and if tu x = tt , m(x).r(x) < . (9) 
where to x is the angle of the boundary at point x. 

Boundary control problem 

Our problem consists in finding To such that for each T > Tq and for every (uq, ui) £ L 2 (fl) x H _1 (f2), 
there exists v £ L 2 (dilN x (0, T)) in such a way that the solution of the wave equation 

' u" - Au = in n x (0, T) , 
u = on dVt D x (0, T) , 

(£) < u = v on dn N x (0,T) , 

u(0) — u in fl , 
•u (0) = u% in SI . 

satisfies 

u(T) = u'(T) = in fl. (10) 
Theorem 3 Assume that (0), ([3]) and (g]) fto/d. 

TTien i/T > 2^%, for every initial data (uo,ui) £ L 2 (J1) x H 1 (£V), i/iere exists a control function 
c(m) 

w G L 2 (3f2jv x (0, T)) smc/i £/ia£ £/ie corresponding solution of (S) satisfies final condition (10). 
Our paper is organized as follows. 

In Section 2, we extend Rellich relations (Theorems [5] and E]) for elliptic problems with mixed boundary 
conditions. 

In Section 3, we apply these relations to prove some stabilization results with linear or quasi-linear 
Neumann feedback (Theorems [l] and 0). 

In Section 4, wc extend some observability and controllability results for the wave equation (Proposition 
Eland ThcoremE]). 

In Section 5, we detail affine examples in the case of a square domain. 

2 Rellich relation 

Here, we briefly extend Rellich relation obtained in [3], [5] to our framework. 
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2.1 A regular case 

We can easily build a Rcllich relation corresponding to the above vector-field m when considered functions 
are smooth enough. 

Proposition 4 Assume that ft is a open set of R" with boundary of class C 2 in the sense of Necas. If 
u belongs to H 2 (fi) then 

2 / Aum..Vudx = / (div(m)l - 2(Vm) s )(Vii, Vit) dx + / {2d v u m.Vu - m.u | Vu| 2 ) da . 
Jn Jn Jan 

Proof. Using Green-Riemann identity we get 

/ Aum.Vu<ix = / d u um.Vuda~ / Vu.V(m.Vu) dx . 

Jn e Jan e Jn E 

So, observing that Vu.V(m.Vu) = im.V(|Vu| 2 ) + Vu.(Vm)Vu = im.V(|Vu| 2 ) + (Vm) s (Vu, Vu), for 
smooth functions u, we get 

2 / Aum.Vucix= / 2<9„u m. Vu da - 2 / (Vm) s (Vu,V«)c(x- / m.V(| Vu| 2 ) dx . 
Jn E Jan e Jn e Jn c 

With another use of Green-Riemann formula, we obtain the required formula thanks to a classical ap- 
proximation. ■ 

We will now try to extend this result to the case of a less regular element u when fl is smooth enough. 

2.2 Bi-dimensional case 

We begin by the plane case: it is the simplest case from the point of view of singularity theory, and its 
understanding dates from Shamir |18j . 

Theorem 5 Assume n = 2. Under the geometrical conditions ([5]) and let u G H 1 (17) such that 

Au e L 2 (Vi) , u /d n D e R 3/2 (dn D ) , d u u /9QN e R 1/2 (dn N ) . 
Then 2d l/ u(m.'S7u) — (m.i/)|Vu| 2 G L 1 (3f2) and there exist some coefficients (c x ) xS r such that 

2 / Aum.Vudx = / (div(m)/ - 2(Vm) s )(V«, Vu) ebc + / {2d v u m.Vu - m.u \ Vu| 2 ) da 
Jn Jn Jan 

+ \ 4(m.r)(x). 
Proof. We follow the proof of Theorem 4 in [5] to get this result. ■ 

Remark 4 As in Theorem ^ of |3J/ ; the assumption 7i (dfln) > is not necessary in the above proof. 

2.3 General case 

We now state the result in higher dimension. 

Theorem 6 Assume n > 3. Under geometrical conditions ([T]) and (J3]), let u G H 1 (f2) such that 

Au g L 2 (Q) , u /ann G n 3/2 (dn D ) , d v u /anN G R 1/2 (dfi N ) . 

Then 2d v u(m.Vu) - (m.v)\Vu\ 2 G L x (9fi) and there exists ( G U 1/2 (T) such that 

2 / Aum.Vudx = / (div(m)i - 2(Vm) s )(Vu, Vu) dx + / (2d v u m.Vu - m.u | Vu| 2 ) da 
Jn Jn Jan 

+ J m.r|C| 2 d7- 

Proof. We exactly follow the proof of Theorem 5 in [5] to get this result. ■ 



3 Linear and quasi-linear stabilization 

Wc begin by writing the following consequence of results of Section 2. 

Corollary 7 Assume that t i— > (u(t), u'(t)) is a strong solution of (S) and that geometrical additional 
assumption ((JJ) if n > 3 (or ([1]) if n = 2) holds, then, for every time t, u(t) satisfies 

2 / Aum.Vwdx< / (div(m)/ - 2(Vm) s )(V«, Vu) cZx + / (2d„u m.Vu - m.iv | Vu| 2 ) da . 
Jn Jn Jan 

Proof. Indeed, under theses hypotheses, for each time t, (u(t),u'(t)) 6 D(W) so that u(t) satisfies 
hypotheses of Theorems [5] or [BJ The result follows immediately from (0) or @ . ■ 

The main tool in the proof of Theorems [TJ [2] is the following result (see proof in [TU]) which will be 

applied with a = — . 

Proposition 8 Let E : R+ — > R + be a non-increasing function such that there exist a > and C > 
which fulfill 

/oo 
£ Q+1 (s)ds < CE(t). 

Then, setting T = CE a (Q), one gets 

ifa = 0, V<>T, <£(0)expfl-| 



ifa>0, Vi>T, E(t)<E(0) 



T + at 



As usual in this context, we will perform the multiplier method to a suitable m. 
Putting Mit = 2m. Vu + au with a a constant to be defined later, we prove the following result. 

Lemma 9 For any < S < T < oo, the following inequality holds 

[ E^ [ ^(div(m)-a)(u') 2 + ((a-div(m))/ + 2(Vm) s )(Vu,Vu))rfxdt 

J S J- S7 



< 



E^ I u'Mudx 
n 



T , P-i F^,* 



E^E 1 / u' Mudxdt 
s Jn 



[ E^ [ m./y((u') 2 - |Vm| 2 - g{u')Mu) da dt . 
J s Jan N 



Proof. We here follow [BJ. 

We Use the fact that u is solution of (S) and we observe that u"Mu = (u'Mu)' — u'Mu' . Then an 
integration by parts gives 



= 



f E^ 1 [ (u" — Au) Mudxdt 
Js Jn 

= E 2 ^- [ u'MudxV [ E^E' [ u' Mudxdt- [ E 1 ^ 1 [ (u'Mu' + AuMu) dxdt , 

L Jn is 2 J s J n J s J n 

Corollary [7] now gives 

/ AuMudx<a / Auudx + / (div(m)/ — 2(Vm) s )(Vu, Vit) g?x + / (2<9„it m.Vii — m.v 
Jn Jn Jn Jan 

Consequently, Grecn-Riemann formula leads to 

/ AuMudx< / ((div(m) - o)I - 2(Vm) s )(Vit, Vu)dx + / (d v u Mu - m.v | Vu| 2 ) da . 
Jn Jn Jan 
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Using boundary conditions and the fact that Vu = d v uu on dflo, we then get 

/ AuMudx.< / ((div(m) - a)I - 2(Vm) s )(Vu, Vu) dx - / m.u (g(u') Mu + | Vu| 2 ) da . 
Jn Jn Jon N 

On the other hand, another use of Green formula gives us 

/ u'Mu'dx= / (a - div(m))(M') 2 + / m.u \u\ 2 da . 
Jn Jn Jan N 

We complete the proof by summing up above estimates. ■ 
Let us now prove Theorems [T] and [2l 

Proof. Following [TU] and [BJ, we will prove the estimates for (uq,Ui) € D(W) which will be sufficient 
thanks to a density argument. 

Using Lemma we have to find a such that div(m) — a and (a — div(m))/ + 2(Vm) s are uniformly 
minorized on J7, that is, almost everywhere on 



(11) 



div(m) — a > c , 
2A m + (a — div(m)) > c , 

for some positive constant c. The latter condition is then equivalent to find a which fulfills 

cssinf (div(m)) > a > ess sup(div(m) — 2A m ) , 

and its existence is now garanted by ([2|). Moreover, it is straightforward to see that the greatest value of 
c such that (fTTj) holds is 

- fess inf (div(m)) — ess sup(div(m) — 2\ m ^\ — c(m) , 
2 V o o / 



and obtained for a = ao := -fessinf(div(m)) + esssup(di 

2\ n p 

With this value ao , we apply Lemma and get 



iv(m) — 2A r 



2c(m) / E~ dt < 



E 1 ^ 1 I u'Mudx 



T p-l 



E^E' u'Mudxdt 



p-i 
E^~ 



Young and Poincare inequality gives 



It follows then 



i'M u dx 



u ((it') 2 - |Vit| 2 - g(u')Mu) da dt . 



< CE(t) . 



2c(m) / E^ dt < C(E z r(T) + E'^r{S)) + C E^ E' dt 
Js J s 

+ { E^ [ m.u((u') 2 - |Vu| 2 - g(u')Mu) da dt . 

J S J dn^ 

Let da m = m.u da. If we observe that E'(t) = — g(u')u' da m < 0, we get, for a constant C > 
independent of E(Q) if p = 1, 



2c(m) f E^ dt < CE{S) + [ E^ [ {{uf - |Vu| 2 - g(u)Mu) da : 
Js Js Jan N 



dt. 
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Using the definition of Mu and Young inequality, we get for any £o > 

2c(m) f E 2 ^ dt < CE{S) + f E 2 ^ [ {{u') 2 (l + ||m||^) + ^-g{u') 2 + e u 2 ) da m dt . 
Js Js Jan N 4e 

Now, using Poincare inequality, we can choose Eq > such that 



/ u 2 da m <-^- [ \Vu\ 2 dx<c(m)E. 



So we conclude 

c(m) f E 2 ^ 1 dt<CE(S)+C [ E 2 ^ [ {{u'f + g(u'f) da m dt . 
Js Js Jon N 

We split dflN to bound the last term of this estimate 

on^ = {x g dn N I |u'(x)| > 1} , dn% = {x e dn N / |«'(x)| < 1} . 

Using ((5|) and ([6|), we get 

f E 2 ^ [ {{u'f + g{u') 2 ) da m dt < C [ E 2 ^ [ u'g(u') da m dt < CE(S) , 
Js Jdn l N Js Jon N 

where C depends on E(Q) if p = 1. 

On the other hand, using (O, ([6|), Jensen inequality and boundedness of m, one successively obtains 

/ {(u 1 ) 2 + g{u') 2 ) da m <C [ ( U ' 5 ( W ')) 2/(P+1) da m <c( f u'g(u') da m ) * < C{-E')^ . 
Jon 2 N Jdn% K Jdn 2 N ' 

Hence, using Young inequality again, we get for every e > 

E*r {{u'f + g{u'f) da m dt < {eE 2 ^ - C{e)E') dt < e E 2 ^ 1 dt + C(s)E(S) . 
s Jdn 2 N Js Js 

Finally we get, for some C{e) and C independent of E{0) if p = 1 

c(m) / E 2 ^ 1 dt < C{e)E(S) + eC E 2 ^ dt . 
Js Js 

c(m) 

Choosing now eC < — - — , one obtains 

c(m) / E 2 ^ 1 dt < CE(S) , 
Js 

and Theorems can be deduced from Lemma HI ■ 

Remark 5 As stated before, we can replace m by Am for any positive A. One can wonder what happens 
to the speed of stabilization 9 = c '™' found in Theorem 3. In fact, a careful estimation of all terms shows 
that one can obtain 

a 2 

C* = fc_ + fc+A 2 + k+-^{\ + C P )C Tr \ 3 , 

where Cp denotes the Poincare constant and Ctt the norm of the trace application Tr : H 1 (O) — * L 2 (<9f2). 
The speed found in our proof is consequently 



c(m) (j- + fc+A + fc+^(l + C P )C T A^j 



It can be shown that 9 reaches a maximum at some point 

fc_ \V3 2 \ /ft. 

A e 



.k + al{l + Cp)C Tr ) 1 a^{l + C P )C T r^ \l k + 
Besides, 9 tends to when A — > or oo. 



Remark 6 In fact, one can replace the feedback law m.u g(u') by a more general one g(yi,u') provided 
that, for some constant c > 1 , 

c _1 (m.i/)p |s| 5+ p < |g(x, s) < c(m.u) p \s\ 5 + p for a.e. x G 8£In and \s\ ^ 1 , 
c _1 (m.i/)|s| < |g(x, s)| < c(m.i')|s| for a.e. x G dQiv and \s\ ^ 1 . 

The details are left to the reader but the previous proof works also in this case. 



4 Observability and controllability results 

It is "well-known that micro-local techniques [3J characterize all partitions of the boundary such that this 
result holds, but constants are not explicit. Thus, using this new choice of multiplier, we will enlarge the 
set of geometric examples with explicit knowledge of constants. We here follow [16] . 



4.1 Preliminary settings 

Following HUM method [13], controlabillity of problem (E) is equivalent to observability of its adjoint 
problem, the solution of the control problem is equivalent to studying the observability properties of the 
adjoint problem. For each pair of initial conditions ((fo, tp\) G Hg(S!) x L 2 (f2), let us consider the solution 
tp of the following wave problem, 



p" - Ap = in il x (0, T) , 

p = ondfix(0,T), 

tp(0) =(p in , 

<p'(0) = p\ in f2 . 



Observability of (E') is equivalent to the existence of a constant C < oo independent of (ipo,ipi) such 
that 



Ea = l f (|^i| 2 + |V^ | 2 )rfx^C / \dM 2 dadt. 

Z JQ Jdn N x(0,T) 

Let us define the operator Wo on Hj(r2) x L 2 (fi) by 



Wo(v,^) = (-^-A V ), 

D(W ) = {{<p, </>) G Hj(fi) x Hj(fi) / Ap G L 2 (ft)} ; 



so that (E') can be written as follows, 



(tp, ip)(0) = (<po,<Pi) ■ 



Remark 7 If (tp,ip) G Z?(Wo), ip is the solution of some Dirichlet Laplace problem and hence regular 
(that is p G H 2 (fi)j. 

Wo is a maximal-monotone operator on Ho(fi) x L 2 (f2) and D(Wo) is dense in Hj(fi) x L 2 (f2) for the 
usual norm. Using Hillc-Yosida Theorem, it generates a unitary semi-group on Hj(Jl) x L 2 (0), we denote 
its value applied at (tpo, p>\) at time t by (tp(t) , ip' (t)) . 
As a consequence, we get conservation of energy. 

Proposition 10 Ift^O and ip is a weak solution of (E')> then 

E(t) = ±J^'(t)\ 2 + \\7p(t)\ 2 )d X = E . 
A weak solution of (E') hence belongs to C(R + ; H^(Q)) n C^R+j L 2 (fi)). 

A solution with (ipo, pi) G £>(Wo) is called a strong solution and satisfies (ip, pi') G C(M+; D(Wq)). 
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4.2 Inverse inequality and exact controllability 

Wc keep similar notations as in Section 2: ao = — (ess inf (div(m)) + ess sup(div(m) — 2A n 

Proposition 11 If T > 2 for each weak solution tp of (£'), the following inequality holds 

c(m) 

ess sup|m.i/| 

E " < ^7 ^t W 9II ll S / \dM 2 d*dt. 
2(c(m)T-2||m|| 00 ) Jan N x(o,T) 

Remark 8 In the case m(x) = (dl + A)(x — xq) with A skew- symmetric matrix, we recover classical 
results (see fTJ ^ .flb J ). 

Proof. Let (ipo, ipi) G D(Wo). We use again M<p = 2m. + aotp. Using the fact that ip is solution of 
(£') and observing that ip"M(p = (cp'Mip)' - ip'Mtp', we get 

f T f r f i T f' T f 

= / / (-cp" + Atp)Mtpdxdt = - / p'Mipdx + / I {p Mp + ApMtp) dxdt . 
Jo Jn k/fi Jo Jo 7n 

As well as in the proof of Theorems [T] and one uses Green- Ricmann formula and Proposition [J] to get 

/ ApMpdx = / ((div(m) — ao)I — 2(Vm) s )(V</?, Vcp) fix + / (d v p>M <p — nx.v \Wp\ 2 )da. 
Jn Jn Jan 

Dirichlet boundary conditions lead to 

/ AipMtpdx = / ((div(m) - a )I - 2(Vm) s )(V^, V» dx + / m.u \d„p\ 2 da . 
Jn Jn Jan N 

On the other hand, another use of Green formula gives us 



/ (p'M<p'dx= j (ao - div(m))|y/| 2 dx, 
Jo Jn 



m Jn 

so, we finally get, using the same minoration as in proof of Theorems Q] and [2] 

c(m) [ [ \<p'\ 2 + \S7p\ 2 dxdt ^ -[ / y/A/^s dx| T + / m.i/ Id^l 2 dcr . (12) 
Jo A2 k/n Jo Jon N 

Using the conservation of the energy, the left hand side in (|12p is 2cTEq. It only remains to estimate the 

r f ~\ T 
term — / p'Mp dx to end the proof. 

L Jn _ J o 

Let us fix a time i G {0,T}. Cauchy-Schwarz inequality leads to 

-J</M<pdx^ (|j^| 2 ) 1/2 (|jil^| 2 ) 1/2 

Denoting by ||.|| the L 2 (fi) -norm, we get the following splitting 

\\Mp\\ 2 = ||2m.v>|| 2 + a 2 1 ||v?|| 2 + 4a / pm.Vpdx. 

Jn 

Green-Riemann formula and Dirichlet boundary conditions give 

j tpm.'Vipdx = — / div(m)|<p| 2 dx, 
Jn 2 Jq 

and since ao — 2div(m) ^ ao — div(m) ^ — c(m) a.e., we finally get that ||M</?|| ^ 2||m|| 00 ||V</?||. 
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Consequently, with Young inequality, we get the following estimate 

- / V )'M ¥ )dx<2||m|| 00 ( / |^f) (/ |V(^| 2 ) ^2\\m\\ oo E . 

So ([12]) becomes 

2(c(m)T - 2||m|| oo ) J B < / m.i/ |<%</f da , 
which ends the proof of Proposition !!!! using the density of the domain. ■ 

Now we can deduce our exact controllability result (Theorem [3]) from Proposition [Tl] following HUM 
method (see [13], Chapter IV). 



5 Example 

Let us consider here the case of a square domain = (0, l) 2 with the following affinc multiplier, 



m(x) 



cot I 



cot 6*9 



(x - x ) 



(13) 



where 6\ and 9 2 belong to ^0, . 

We will discuss the dependence of d^N and dilo on x . 

First let us consider one edge [ab] of with its normal unit vector i>. One can easily see that 

1 



m(x).i/(x) 



sm 



-(X-X O ).U0. 



where = 9\ (resp. 6 2 ) if [ab] C [0, 1] x {0, 1} (resp. {0, 1} x [0, 1]) and vg is deduced from v by rotation 

of angle —6. Without any restriction, we suppose a.vg < b.vg. 

Then there exists an interface point along [ab] if and only if Xo belongs to the belt 

Bg = {x e R 2 I EL.Ug < X.Ug < b.Vg} . 

In this case, at this interface point Xi, we get with similar notations, 

m(xi).r(xi) = — -(xi - x ).x 9 . 

smu 

Then additional geometric assumption ([7]) is not satisfied if xo belongs to half-belt B^ (see Fig. 1). 

B- 




Figurc 1: If xo belongs to B g , we get mixed boundary conditions along [ab] and condition ([7]) is satisfied. 
We now can describe every situation by considering only three following cases (see Fig. 2), 



(CI): 0<6 1 <6 2 <^ 1 (C2): O<0i<|<0 2 <| 



(C3): ~<0i<0 2 <|. 



We also show a fully detained partition in some particular case coresponding to (C2) (see Fig. 3). 
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Figure 3: Example of Dirichlet and Neumann parts of the boundary in a case of (C2)-type. 
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